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Abstract 

We present a generic Z2 x Z2-invariant scalar field theory with four real scalar fields in six-dimensional Minkowskian spacetime 
which yields solutions consisting of two intersecting domain- wall kinks which are each paired by fields with lump-like profiles. 
For a special parameter choice, analytic solutions can be obtained. We show that the Z2 x Z2 symmetry can be maintained 
while coupling fermions by introducing scalar Yukawa couplings to one kink-lump pair and six-dimensional pseudoscalar Yukawa 
couplings to the other, and we show that there exists a zero mode localized to the domain-wall junction in this case. We also 
show that scalar fields can be localized. 



PACS numbers: 

I. INTRODUCTION 

The possibility that there could be hidden extra dimen- 
sions of space continues to be an intriguing one for de- 
scribing physics beyond the Standard Model (SM). Most 
extra-dimensional models, such as string theory, use com- 
pact and tiny Planck-scale extra dimensions. Recently, 
large extra dimensions of order an inverse TeV have been 
considered ^-E] , which became particularly popular with 
the advent of the model by Arkani-Hamad, Dimopolous 
and Dvali (ADD) which featured a solution to the hier- 
archy problem. Later, Randall and Sundrum proposed 
the RSI model 7 which presents an alternative solution 
to the hierarchy problem in a slice of warped anti-de Sit- 
ter (AdS) geometry. The same two authors showed that 
infinite extra dimensions were possible, in a model now 
termed RS2 [Sj. In this model, matter is confined to a 
fundamental positive-tension brane put into the theory 
by hand, and gravity is localized onto this fundamental 
brane. 

One could consider the possibility that such a brane 
on which matter and gravity is localized is not put in 
by hand but rather generated dynamically by the under- 
lying theory. Such an object must be a stable solution 
to the theory; thus the obvious candidates for such an 
object are solitons. Topological solitons are classical so- 
lutions to the dynamical equations of a theory which non- 
trivially map a boundary of one or more dimensions of 
spacetime at infinity to the moduli space of vacua. These 
mappings belong to the non-trivial homotopy classes of 
mappings between these two spaces. Such topological 
solitons are stable against decay to a trivial vacuum due 
to the existence of conserved topological charges which 
are associated with these homotopy classes. Examples 
of solitons are domain walls which are classical solutions 
which map spatial infinity along both directions of a par- 
ticular dimension to two distinct but degenerate vacua, 
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strings which are non-trivial maps from the boundary of 
a two-dimensional space to a moduli space with the ge- 
ometry of a circle, and monopoles which are non-trivial 
mappings from a sphere at spatial infinity to a sphere 
of vacua. Localizing matter on solitons, particularly do- 
main walls, is not a new idea and the proposition of the 
world volume of a domain wall forming our observable 
3+1-dimensional (3+1D) universe was first proposed by 
Rubakov and Shaposhnikov [jj]. 

In addition to the original proposal by Rubakov and 
Shaposhnikov, there have been several other attempts 
at demostrating the viability of such a model based 
on domain-wall branes. Localization of fermions and 
scalars can be easily achieved by introducing appropri- 
ate Yukawa and quartic scalar interactions respectively 
(for an analysis of the localization of scalar and fermionic 
zero modes as well as their Kaluza-Klein (KK) spectra, 
see [10]). For the localization of gauge bosons, Dvali and 
Shifman have proposed a mechanism by which a gauge 
group G is broken down to a subgroup H in the interior 
of the domain wall and is unbroken and confining in the 
bulk, localizing the gauge bosons of H to the domain wall 
through confinement dynamics |llj . Importantly, gravity 
can be localized on the domain wall, yielding an RS2-like 
warped background geometry and in its presence, scalars 
and fermions can still be localized [12]. A viable model 
utilizing the Dvali-Shifman mechanism with G = SU(5) 
and H = SU(S) x SU(2) x U(l) in 4+1D spacetime was 
proposed in [13] which has some very interesting phe- 
nomenology; the SM components embedded within the 
SU(5) multiplets are split and localized about different 
points in the extra dimension leading to a natural re- 
alization of the split fermion mechanism first proposed 
by Arkani-Hamed and Schmaltz [14] . Furthermore, such 
splitting of fermions and scalars within the model can ac- 
count for the fermion mass hierarchy, quark mixing and 
the suppression of proton decay [15j and a simple exten- 
sion of the model to include the discrete flavor group A4 
can naturally generate large lepton mixing angles [16] . 
Models based on extended gauge groups such as 50(10) 
[17] and E 6 [18] are also possible. In this paper, we pro- 
pose an extension involving the addition of a second extra 



dimension with matter being localized to the intersection 
of two domain walls. 

Topological defects of co-dimension 2 have been con- 
sidered as candidates for the localization of helds onto a 
3+1D subspace. The simplest example of a co-dimension 
2 defect is a string and it has already been shown that 
3+1D gravity can be reproduced on such an object in a 
5+1D spacetime [THIHO]- Extensions of the RS2 model 
for the intersection of n fundamental branes in 4+n- 
dimensional spacetime were proposed in |21l I22j . Other 
than strings, we could instead consider the possibility of 
using a pair of domain walls in 6D spacetime to localize 
fields on to a 4D world volume. 

There are two ways of introducing a second domain 
wall in order to freeze out a second extra dimension. One 
way is to localize a scalar field onto a domain wall and 
have that scalar field develop a tachyonic mass which in- 
duces a breaking of a second discrete symmetry, yielding 
a second domain wall localized to the first. This is called 
a nested domain wall or domain ribbon and there is some 
literature that has dealt with this scenario [231 HI] an d 
with the possible localization of gravity on such a defect 
[25j . A second way to produce a dimensional reduction 
from 5+ ID to 3+ ID is to have two stable domain walls 
which intersect. There have been some supersymmctric 
models which yield rare exact solutions for a pair of non- 
trivially intersecting walls |26j . For more on models with 
intersecting domain walls or domain-wall junctions, see 



In this paper, we present a Z2 x Z2-invariant 5+1D 
interacting scalar field theory with four real scalar fields, 
where the 6D masses are tachyonic for two of the scalar 
fields and positive definite for the other two, in which a 
rare analytic solution for two intersecting domain walls 
can be obtained for a particular parameter choice. We 
find that there exists a class of energy degenerate solu- 
tions with two domain walls: one in which the walls are 
perpendicular, a range of solutions where the walls inter- 
sect at an angle between and 90 degrees and another 
in which the walls are parallel. We give topological argu- 
ments for why the perpendicular solution cannot evolve 
to the parallel solution given for the particular parameter 
choice, as well as an argument for why the perpendicular 
solution might be energetically favorable to the solutions 
with intersection angle less than 90 degrees in a nearby 
region of parameter space (assuming they are not topo- 
logically distinct). We also show that chiral fermions and 
scalars can be localized to the intersection of the domain 
wall. This is important because these chiral fermions 
and scalars form the building blocks of the quarks, lep- 
tons and Higgs bosons of an effective Standard Model-like 
theory dynamically localized to the intersection. 

In Sec. [IIJ we give a review of domain walls with an ad- 
ditional scalar field which attains a tachyonic mass in the 
interior of the domain wall and which thus condenses to 
attain a lump-like vacuum expectation value profile. In 
Sec. |III| we outline the model generating the intersecting 
domain wall solution, give the form the solution for which 



the walls are perpendicular and we also give some topo- 
logical arguments for why this solution cannot evolve to 
solutions for which the walls are either intersecting at an 
angle between and 90 degrees or parallel despite these 
solutions being energy degenerate. In Sec. IV we discuss 
fermion localization and we show that localization of a 
single chiral zero mode on the intersection of the domain 
walls is possible. In Sec. |Vj we show that scalars can also 
be localized to the domain- wall intersection. Section IVTl 
is our conclusion. 



II. DOMAIN WALLS WITH INTERNAL STRUC- 
TURE: A REVIEW 



Before we do the full analysis for the intersecting 
domain-wall solution, we present a brief overview of 
domain-wall solutions with internal structure. Domain 
walls with internal structure are common in interact- 
ing scalar field theories involving at least two real scalar 
fields. They consist of one scalar field which generates 
the domain wall and thus has a kink-like profile, and one 
or more scalar fields which interact with the kink field at- 
taining a lump-like profile. One can think of these extra 
fields as being dynamically localized to the wall which 
at the same time attain tachyonic masses in the interior 
of the wall and thus non-zero vacuum expectation values. 
In general, if the back reactions of these fields on the kink 
field are significant, these can affect the domain- wall so- 
lution, particularly its width. 

To present an example of a domain wall with internal 
structure, we write down a Z2-symmetric scalar field the- 
ory with two scalar fields, 77 and x- The field rj will form 
the background domain wall, and x will condense in the 
interior. We desire to find a solution such that rj inter- 
polates between two values ±tj from y = —00 to y = +00 
and where the vacuum expection value of x tends to zero 
out at infinity. To generate such a solution and ensure 
its topological stability, we must ensure that 7/ = ±v, 
X = are global minima of the potential. Hence, we give 
77 a tachyonic mass, while the mass squared of x will be 
made positive definite. Given these requirements and the 
Z2 symmetry 77 — > —rj, x ™>> ~Xi the scalar potential is 



V(V,X) = \\{V 2 v 2 ) 2 + \\M - v*)x 2 + \4x 2 
+ \ X xX A + 9 V x r fx + KxVX 3 - 

(1) 

To ensure that the potential is bounded from below as 
well as to make sure we have the desired properties for 77 
and x, we make the following parameter choices, 

A,, > 0, A x > 0, X, lx v 2 > p? x . (2) 

To form the kink-lump background, we need to look 



2 



for a static solution to the Euler-Lagrange equations 



□?/ + A r) (?7 -v )ri + \ 11x rix + ig nx ri x + KxX = 0, 
Ox + A x x 3 + M X X + x vx(v 2 ~ v 2 )x + g n xV 3 + 3h vx T] 2 x 



= o, 



subject to the boundary conditions 
v(y = ±°°) = ±f) 

x(y = ±°°) = o, 



(3) 



(4) 



that depends solely on the coordinate y so that r\ = rj(y) 
and x = x(y)- This equation can be solved numerically 
and in general the profile for ry is kink-like and the pro- 
file for x is indeed lump-like. For the special parameter 
choice 

9vx = Kx = °: 2 Mx(V- A x) + (\ A x- A to) w2 = °' ( 5 ) 
one finds the analytic solution 



v(y) — v tanh (ky), 
X(y) = Asech (ky), 



(6) 



where k 2 = fj, 2 and A 2 = XnxV 2m * . A plot of the solu- 



tion is given in Fig. [T] 
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FIG. 1: A plot of the profiles for 77 and x- 



III. THE BACKGROUND SCALAR FIELD THE- 
ORY 

Our model describing the background fields is a 5+1- 
dimensional (5+ ID) quartic scalar field theory with four 
scalar fields invariant under a Z2 x Z2 symmetry. Under 
this symmetry we assign the following parities to these 
fields: 



m 



(-.+) 
(+,-) 



xi ~ (-,+), 
X2 ~ (+,-)• 



(7) 



The fields r\\ and 772 will form the two perpendicular back- 
ground domain-wall kinks while xi and X2 will attain 
lump-like profiles parallel to each respective wall. 

Given the parity assignments, we may write the most 
general quartic scalar potential of this theory as 



Vdw = (rjl - vl) 2 + -X mxi (17? - vl)xl + ^xiXi 
+ \ X xiXl + QmxiViXi + h mxl ViX 3 i 

. 1 x / 2 2\2 ,1% f 2 2\ 2 , 1 2 2 

+ - Vfe - v 2) + v*V2X2 (% - V 2 )X2 + ~M X2 X 2 



+ 1K2X2 + 9mx2 1 l 2 X2 + h V2X2 mX 2 



+ \ X ximXi(vl - v\) + \\ lX2 x\x\ 

+ oVxi^mXi^l + 2 A ')ixiX2'7iXiX2 
+ A »nxi^2X2 7 7iXi ? 72X2- 

(8) 

Choosing parameters such that this potential is un- 
bound from below, including X Vl , \ m , A Xl , A X2 > 0, there 
are four global minima given by r\\ = ±«i, r\2 = iu2j 
Xi = X2 = 0. 

To set up an intersecting domain wall solution, we must 
find a static solution for all of the four fields rji, 772, %\ 
and X2 to the Euler-Lagrange equations which at the very 
least interpolate amongst the four vacua at infinity along 
the corners of each quadrant in the y — z plane. Firstly, 
we attempt to find a solution in which the walls are mu- 
tually perpendicular. Along one edge where 771 is con- 
stant at one of its vacua ±«i and where also xi = 0, the 
field 772 should interpolate between the values ±tj 2 and X2 
should be zero at infinity along the edge and condense in 
the middle of the edge, much like the one-dimensional 
domain wall with internal structure discussed in the pre- 
vious section. The same should apply to 771 and xi along 
the edges where 772 and %2 are fixed. This motivates us 
to look for solutions obeying boundary conditions of the 
type 

771(7/ = ±00, z) = ±vi, 771(7/, z = ±00) = vi tanh (ky), 
772(2/ = ±00, z) — v 2 tanh (Iz), rj 2 (y, z = ±00) = ±tj 2 , 
Xi(y = ±00, z) = 0, xi(y, z = ±00) = Ax sech (ky), 
X2(y = ±00, z) = A 2 sech (Iz), xa(y, z = ±00) = 0. 



By making the parameter choice, 



A m»?2X2 — A Xl>)2X2 — A »ilXl>j2 — A >JlXlX2 — A »ilXl>j2X2 

9mxi = h"nixi = .9»?2X2 = ^>j 2 x2 = 0) 
^vim v i = X XlV2 A 1 , A r;iI?2 f 2 = X mX2 A 2 , 

A »7lX2 U l — A XlX2^1i A Xl»72 U 2 = A XlX2^2; 



(9) 



0. 



2 /^ Xl ( A »nxi ^ A xi) + ( A m A xi _ ^ixi)" 2 - °' 
2/-t X2 (A, )2X2 - A X2 ) + (A^ 2 A X2 - A^ 2X2 )u 2 = 0, 



(10) 
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one finds that 

Vi(y) = v i tanh(%), 
rj 2 (z) = v 2 ta,nh (lz), 



Xi(y) = A 1 sech(ky), 
Xi{z) = M sech (lz), 



(11) 



where k 2 



Xl 



I 2 



-2M 



is a solution to the four coupled Euler- 
Lagrange equations resulting from the potential in Eq. 
[8] and satisfies the boundary conditions in Eq. [9j We 
give plots for r/i, Xii 772 and \2 for the solution in Eq. |11| 
in Fig.[2j Fig.[3j Fig.|4]and Fig. [^respectively, in terms of 
the non-dimensionalized coordinates y — ky and z = lz. 
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FIG. 2: A plot of 771 for the solution in Eq. 11 
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FIG. 4: A plot of 772 for the solution in Eq. 11 
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FIG. 5: A plot of X2 for the solution in Eq. [TTj 




FIG. 3: A plot of \i f° r the solution in Eq. 
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It is important to note that there are other solutions 
with two kink-lump pairs to the coupled Euler-Lagrange 



equations. Without loss of generality, assume that 771 and 
Xi are the same as in Eq. [TT] but consider instead that 
772 and X2 take the form 



772(7/, z) = i'2 tanh [l(cos9y + sin#z)], 
X2(y, z ) = A 2 sech [I (cos By + sin 9z)\. 



(12) 



Acting with the dAlembertian operator on 772 and X2 
yields 



2 (vi - 



2l 2 

□772 = - 

□X2 = ^2 Xl - l2 X2- 

The resulting relations of the kinetic terms of 772 and 



(13) 



X2 to themselves given in Eq. 13 are independent of the 
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relative angle, 9, between this kink-lump pair and that 
for rji and xi- Hence, under the same parameter choice 
as in Eq. (TUJ 



Vi(y) = -Manh(fcy), 

Xi(y) = A i sech(fcy), 
r\i{y,z) — v 2 tanh [I (cos 9y + sin 9z)] , 
%2(y, z) = A 2 sech [Z(cos6»y + sinflz)], 



(14) 



is a solution to the Euler-Lagrange equations for gen- 
eral 9. This means we have a whole class of solutions 
ranging from a solution in which the two walls are paral- 
lel (9 = 0), through intermediate angles of intersection, 
to the perpendicular solution of Eq. [TT] The parallel and 
angled solutions obviously do not satisfy the same bound- 
ary conditions as the perpendicular solution given in Eq. 
[9j The angled solutions still divide the y — z plane into 
four domains which tend to all four of the discrete global 
minima out at infinity. The parallel solution is a wall be- 
tween two of these discrete vacua out at infinity along one 
direction. In calculating the energy density, the kinetic 
terms and the potentials describing the self-interactions 
and mass terms involving only r]i and xi an d likewise 
only r/2 and xi are associated with the energy density of 
the single kink-lump pairs, and integrating them over the 
directions normal to the walls yields the tensions associ- 
ated with each kink-lump pair. In general, in theories 
with multiple domain-walls, the quartic interactions be- 
tween the fields yielding different domain walls lead to 
a tension associated with the intersection or junction of 
those walls. In this scalar field theory and with the pa- 
rameter choices we made in Eq. |10| we find that this 
tension is precisely zero, 

dm = ^-Mslli - vl){ril - vl) + ^X mX2 (vl - vl)xl 

|1\ 2/2 2\,1\ 22 

+ 2 A xi%Xi(% - v 2) + 2 A xiX2XiX 2 > 
= - \ 1 i- l2 v i v 2 sech 2 (ky) sech 2 (lu) 



1 

2 
1 

2 
1 

2 
= 



\iX2 v \A\ sech 2 (ky) sech 2 (lu) 
X Xim A\v\ sech 2 (ky) sech 2 (lu) 



+ ^ A xiX2^i^2 sech 2 (ky) sech 2 (lu), 



(15) 

for all angles 9 where here we have used u = cos 9y + 
sin 9z. These solutions are degenerate in energy. 

It turns out that despite the energy degeneracy, nei- 
ther the perpendicular nor angled solutions can evolve 
into the parallel solution. This is not surprising as the 
former two interpolate amongst the four vacua along the 
boundary, while the latter only interpolates between two 
of them so we expect them to be topologically distinct. 
To be precise, there exists a topological charge associated 
with the 2-dimensional boundary of the y — z plane. The 



associated topological current is defined by 

jMNOP = e MNOPQR e ij dQrj . dRrj , (lg) 

Clearly, the 6-divergence of this current vanishes and is 
thus conserved. The conserved topological charge asso- 
ciated with this current is just 



Q 



ABC 



= / <f X J' 



6„ tO ABC 



(17) 



Since the background fields are solely dependent on y and 
z, only the elements Q l ^ k where i, j, k = 1, 2, 3 are non- 
zero. One can show that these charges are proportional 
to the integral 

1=1 dydz(d i r] 1 d 5 ri2 - 8 5 r]id4ri2), (18) 

where S denotes the y — z plane. Using Stokes' theorem, 
one can also write this as 



1=1 r]id A r] 2 dy + rjid^dz, 

'9S 



(19) 



77iV?72.dl. 



OS 



One can easily show from Eqs. [T8| and [T9| that for the 
perpendicular and angled solutions I = 4:V\V2 whereas for 
the parallel solution 1 = 0. Thus the perpendicular and 
angled solutions are stable against decay or evolution to 
the parallel solution despite the energy degeneracy. 



The topological charge in Eq. 17 does not differentiate 
between the perpendicular and angled solutions. It is not 
clear whether this means that these solutions are in the 
same topological class as there could exist other topologi- 
cal charges which differ between the two types of solution. 
In case they are topologically equivalent, one can imag- 
ine that one could perform a small perturbation from the 
parameter region considered in generating the analytic 
solutions to ensure that the perpendicular solution is the 
most energetically favorable one since the energy degen- 
eracy between the solutions is likely not true in general. 
From a rough glance at the interactions in Eq.[l5]one can 
see that if one performs a perturbation X mm — > X m rj 2 + e > 
with e > 0, that the contribution e?y 2 ?7| is minimized for 
9 = 90° and tends towards infinity as 9 — > 0. Unfortu- 
nately, there are also resultant perturbations to the fields 
Vii V2, Xi an d X2; and these perturbations satisfy four 
coupled partial differential equations which can only be 
solved numerically, so we can not give a definitive answer 
here. 

We also need to perform a local stability analysis of 
the solutions given in this section. Likewise, this requires 
numerically solving four non-linear coupled partial differ- 
ential equations and we defer this to a later study. For 
the rest of this paper, we will assume that the perpen- 
dicular solution is stable and that one can always choose 
this to be the background solution. 
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IV. FERMION LOCALIZATION 

In this section, we show that fermions can be localized 
to the intersection of the two domain walls. Normally, in 
the case of a single domain wall, one localizes fermions 
to the centre by Yukawa coupling them to the relevant 
scalar field. In the case of a single domain wall with a 
Yukawa coupling of the form ^/^S 77, as r\ — > — r/ under 
the discrete Z2 symmetry, to preserve the symmetry and 
maintain topological stability, we require that the Dirac 
bilinear — > — under the symmetry. This can 
be achieved by choosing the fermionic fields to transform 
individually as VP — ¥ ±iT^, where T is the gamma matrix 
associated with the direction parametrizing the profile of 
the domain-wall solution. In 4+1-dimcnsional (4+1D) 
theory, one often chooses r = T 4 = — i^ 5 for example. 

In six dimensions, with our background set-up, if one 
wishes to localize chiral fermions to the intersection re- 
gion of the background solution, one must Yukawa cou- 
ple the desired fermionic fields to all four scalar fields. 
However, with two independent Z2 symmetries, one has 
a problem in attempting to couple a fermionic field to 
both defects since if we perform the first Z 2 transforma- 
tion rj! -t -Tii, xi -^_rXi, m -> V2, X2 -> X2, and 
maintain that — > — 4"!' under such a transformation, 
then \E"]/T72 and \E"I'X2 are not invariant under the first 
symmetry, and likewise tytyrji and "J^Xi won't be if we 
impose the second Z2 symmetry. Hence, we cannot lo- 
calize fermions to the defect by using just scalar Yukawa 
interactions to all fields without compromising topologi- 
cal stability. This reflects the fact that we must choose 
our effective 6D localization bulk mass matrix carefully 
in order to localize chiral fermions on the intersection, an 
issue that was first raised in |22j . 

However, in 5+lD there is another possibility, since in 
spacetimes of even dimensionality there always exists a 
chirality operator and thus there always exists a pseu- 
doscalar bilinear in these spacetime dimensionalities. In 
5+lD, the chirality operator T 7 is defined 



p7 -pi -p2-p3p4-p5 

1 



6! 



e MN0 PQRr M T N T T p TQT R . 



(20) 



One can then define the pseudoscalar bilinear 'IT 7 'I'. 
Now we consider the Yukawa terms 



C Y uk = -ih m ^T 7 ^r]i - ih Xl ^T 7 ^ X i 



(21) 



and ask if it is possible to define two independent trans- 
formations for each Z2 symmetry for \P such that for the 
first symmetry in which r/i — > —r/i and xi ~~ ► ~Xi we 
have *r 7 * -> - *r 7 * but -> while for the 
second reflection symmetry r]2 — > —772 and X2 — > — X2 
we have ^T 7 ^ unchanged but ¥\I> -> -W*. Due to 
the fact that T 7 anticommutes with the gamma ma- 
trices, for the second Z2 one can easily show that the 



usual type of transformation —> iT 5l $> can be chosen. 
For the first Z2, one can show that the transformation 

-> ir 4 r 7 * induces the transformation _> ^vp and 
*r 7 * -> -*r 7 *. Hence, we have shown that there ex- 
ists a mechanism to couple a fermionic field to all four 
background scalar fields with the combination of scalar 
and pseudoscalar Yukawa couplings given in Eq. |21| 

Henceforth, we assume boundary conditions such that 
we have the perpendicular intersecting domain-wall so- 
lution and we take the following sets of transformations 
to be our reflection symmetries which ensure topological 
stability of the background, 



y-> -y, 

z — > z, 

m - vi, 
xi - Xi, 

V2 -> V2, 
X2 -> X2, 



(22) 



and 



y -> y, 

z — > — z, 
Vi -> Vi, 
Xi Xi, 

m -> - ?n, 

X2 -> - X2, 



(23) 



We now need to show that there is indeed a chiral zero 
mode localized to the intersection of the domain walls. 
Writing down the resultant 6D Dirac equation for <]/, we 
have 



iT M d M ^ + iWi(y)T 7 V- W 2 (z)^ = 0, (24) 



where 



Wi(y) = h m rn(y) + h XlX i{y), 
W 2 (z) = h V2 T] 2 {z) + h X2 X2{z). 



(25) 



In order to perform dimensional reduction and calcu- 
late the profiles of the modes of \P, we must choose a 
basis for the 5+lD Clifford algebra. One can show that 



r 4 = 01 ® - i-f 



-ia 3 <g> 1 



— %x 
vy 5 

-i 0' 
% 



(26) 



satisfies the 5+ ID Clifford algebra 

{T M ,T N } = 2 V MN 1 



(27) 
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and is thus an appropriate choice of basis for the 5+1D 
gamma matrices. In this basis the 6D chirality operator 
is 



CT 2 ® 1 



% 

-i 



(28) 



Decomposing 'F into components ^>± which have 4 com- 
plex components and are eigenvectors of T 5 , 



(29) 



one can shown that the 6D Dirac equation reduces to 
(iVfy + 7 5 a 4 + - 3 5 *- - VF 2 (z)<F_ 



= 0, 
(30a) 

W 2 (z)^+ = 0. 

(30b) 



To calculate the profiles of all modes, due to the fact 
that the excited Kaluza-Klein (KK) modes are usually 
Dirac fermions it is useful to find the corresponding 
Klein-Gordon (KG) equation that the components ^± 
satisfy. Op erati ng with (ij^d^ +7 5 c>4 — W\(y)) from the 



left on Eq. 30a and likewise {i^^d^ + ^ d± + W\{y)) on 



Eq. |30b[ one obtains the KG equations 

n^ ± +{W l ( y ) 2 TW[{y) lb )^ ± +{W 2 {z) 2 ±W^ z ))^ ± = 0. 

(31) 

Now we expand each of *± as a series of modes. As one 



can see from Eq. 30 the y-dependent piece is chirality- 
dependent while the z-dependent piece is not. The z- 
dependent piece is only dependent on whether the com- 
ponent is \I/ _ or Thus we make the expansion 



*±(a; M ,y,z) = 

, ±L (y)g^(z)^M + f? R (y)g n ± l (z)^ R (x,) 



(32) 

Here, m just denotes some mass eigenvalue, ip± L R de- 
notes the 3+1D left/right-chiral mode of mass m embed- 
ded in the component <F± of 'J', and f± L R (y) and g±(z) 
are profiles for these modes along the y and z directions 
respectively. Since any 3+ ID fermionic mode should sat- 
isfy a corresponding Klein-Gordon equation, let us sub- 
stitute the expansion into Eq. |31| and demand that modes 
satisfy n 3+1 ip^ LR (x f _ l ) 
Eq. [31] reduces to 



for the left and right-chiral components respectively. 

Demanding that the profiles satisfy the following 
Schrddinger equations (SE) 



dy 2 

a J±R 

dy 



+ (W^y) 2 ± W[(y))fZ l L (y) = Xij^iy), 

(34a) 

, + (W^yf T W[(y))f¥ R (y) = A^/S^), 

(34b) 



d 2 g% 
' dz 2 



+ {W 2 {z) 2 ±W! 1 {z))g^{z) = \ 2 ± gl{z) 1 



(34c) 



we find that the values of the squared masses of the lo- 
calized KK modes are 



l ±L,R 



— A±L o + A 



2 



(35) 



for the modes embedded in ty±. 

Given the definitions for W\ and W 2 in Eq. [25] one 
can see that the potentials of Eq. [34] are hyperbolic Scarf 
potentials. These potentials are well known and can 
be solved analytically [31 33 . For simplicity, let us as- 
sume that both h Vl and h V2 are positive definite. Non- 
dimensionalizing all variables and parameters except m 
as 



h Vl — 



V = ky, 
h m vi 
k ' 
h m v 2 



= I 



z, 



A ±L.R 



A ±L,R 
k 2 ' 



h Xl = 



h-X2 — 



k ' 

h X2 V2 

I ' 



(36) 



X 2 

\2 _ A ± 



and non-dimensionalizing the profiles as 

~ 9 2(z) = H<#(*), 



(37) 



ff±L 

dy 2 



one can show that each of the Eqs. |34a[|34b and 34c has a 
finite number of localized, square-normalizable solutions 
as well as a delocalized continuum. First, let us start 
171 <P±l,r\ x »)- We nnd tnat with the z-dependent equations. For positive h, l2 , it is 

the potential for W 2 — W' 2 which generates a series 
of \h m ~\ localized modes starting from an eigenvalue of 



(W 1 (y) 2 ±W{(y))f^ L (y))]g^(z) 



0. The eigenvalues of these modes are given as 



d 2 g± 



+ .f± L (y) 

= m 2 f : 



1 dz 2 

-_ L (y)g±(z), 



(W 2 (z) 2 ±W'(z))gl\z)] 



X 2 _ n = 2nh V2 



d 2 f. 



(33) 



±R 



dy 2 

l 2 g™ 
m 2 f^ R (y)g^(z), 



+ f± R (v)\ 



(Wi(y) 2 TWl(y))fZ R (y)]gZ(z) 
(W 2 (z) 2 ±Wi(z))g™(z)} 



■, n = 0,l,...,LVJ- (38) 
The A^q profile is given by 

g° (z) = D° er~ hr >2- log f cosh ( 2 )l _2 ' i x2 arctan I tanh (z/2)] 

and the profiles for the excited localized modes can be 
generated by applying the ladder operator which is pro- 
portional to W 2 (z) — 4s- As for the potential for 
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Wf + W2, given our parameter choice there is no solution 
with A+ = 0. Rather, there are [h V2 \ localized modes 
starting from an eigenvalue — 2h rj2 — 1 (provided 
h„ 2 > 1). All this implies that, in considering just the 
interaction with the z-dependent kink-lump, for h V2 > 0, 
there is a massless 4+1D Dirac zero mode generated in 
\I r _ and then a tower of massive 4+1D Dirac modes em- 
bedded in both and In then further considering 
the interaction with the y-dependent part of the solution 
and calculating the J±l,r, each of these 4+1D modes will 
generate a tower of 3+ ID left- and right-chiral modes lo- 
calized to the domain-wall intersection. Obviously, any 
chiral zero mode produced on the intersection must be 
embedded in the massless 4+1D zero mode of \P_ in this 
case. If we choose h V2 to be negative instead, the roles of 
ty- and in this situation are reversed, and the chiral 
zero mode must therefore be embedded in ^! + . 

Now let us analyse the y-dependent equations. These 
equations are also SE's with the hyperbolic Scarf poten- 
tials. In this case the particular form of the hyperbolic 
Scarf potential for the modes is dependent on chirality. 
Since we are assuming h„ 2 > without loss of generality, 
let us focus on the left- and right-chiral modes embedded 
in first, since as noted above this is the component 
containing any potential chiral zero mode. Looking at 
Eqs. |34a| and |34b[ one sees that the potentials for the left- 
and right-chiral components in are W\(y) 2 — W[(y) 
and Wi(y) 2 +W[(y) respectively. Hence, we easily deduce 
that for h Vl > 0, the equation for the left-chiral modes of 
has the same form as that for the z-dependent profile 
equation for and thus has a mode starting from an 
eigenvalue of A 1 l £0 = 0, and likewise that for the right- 
chiral modes has the same form as that for the g™(z) 
and thus only has solutions with positive definite eigen- 
values. Since from Eq.[35]we know that these eigenvalues 
directly contribute to the mass, this implies that for the 
choice h Vl> Q, /i,j 2 >o, there is a single massless left-chiral 
zero mode embedded in localized to the intersection. 

The resultant solution for the A^q = eigenfunction, 

I-l{v) is g iven b y 

f° l(v) = C° L e~' l ' n lo S I cosh (5)]- 2 ' l xi arctan [tanh (y/2)] 

(40) 

and thus the full profile over the y — z plane for the left 
chiral zero mode is 

F _ L (y,z)^f°- L (y)9-^)- (41) 

and a plot of this is shown in Fig. [6] 

The y-dependent profiles for the higher localized left- 
chiral modes can be accessed by applying a ladder opera- 
tor proportional to W\ (y) — The eigenvalues of these 

/™ L profiles are given as 

\ 1 _ Ln = 2h ril n-n 2 . (42) 

Putting Eqs. [38] and [42] together and converting back 
to dimensionful variables, we can see that the resultant 



/ - L (y)g -(*) 

2 




y 2 



FIG. 6: A plot of the profile for the left chiral zero 
mode embedded in >]/_ for the parameter choice 
h Vl = 10, h Xl — —5, h V2 = 20, and h X2 = 4 




FIG. 7: A plot of the profile for the left-chiral i = 1, 
j = mode in for the parameter choice h rn = 10, 
h Xl = —5, h, l2 = 20, and h X2 = 4 

squared masses of the left-chiral modes embedded in 
are 

m 2 _ Li j = }}_ Li k 2 + ^-jl 2 , 

= lih^vxk - 2i 2 k 2 + 2jh V2 v 2 l ~ 2j 2 l 2 , (43) 
i = 0,1, LVJ, j = 0, 1, [h V2 \- 

For the parameter choice h m = 10, h Xl = —5, h m = 20, 
and h X2 — 4, we give a plot of the full profile F°_ L (y, z) 
for the zero mode (i = 0, j = 0) in Fig. [6] We also give 
plots for the analogous profiles of the left-chiral i = 1, 
j = 0, i = 1, j = 2, and i = 2, j = 2 KK modes in Fig. 
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FIG. 8: A plot of the profile for the left-chiral i — 1, 
j = 2 mode in for the parameter choice h Vl = 10, 
h Xl = —5, h V2 = 20, and h X2 = 4 



f 2 - L Cy)g 2 _Cz) 




FIG. 9: A plot of the profile for the left-chiral i = 2, 
j = 2 mode in for the parameter choice h r)1 = 10, 



-5, hri 2 



20, and 



[7j Fig. [8] and Fig. [9] respectively. 

For the right-chiral modes of ^_ , there is no zero mode, 
but the massive modes have the same mass as the left- 
chiral counterparts, and the masses are 



-Rij 



Ai fli fc 2 + A_^Z , 
= 2ih m vxk - 2i 2 k 2 + 2jh V2 v 2 l - 2j 2 l 2 , (44) 
i = l,...,[h m \, j = 0,1,..., [h n2 \. 



For the left-chiral and right-chiral localized modes em- 
bedded in the eigenvalue associated with the z- 
depedent profiles, A 2 ^-, is always more than zero, thus 
all these modes are massive. Also, A 2 ^- = A 2 ^ for 



j = 1,2,..., J ■ If one looks at Eqs. (34a) and 34b 



one sees that the y-dependent profiles /™ L of the chi- 
ral modes of satisfy the same equation as those of 
equal mass and opposite chirality in \& + , In other words, 
f±L = which in turn implies that \ 1 +L R — X_ R L . 
Thus the masses for these modes are simply 



l +Lij 



+:i L 



= 2ih m vik - 2i 2 k 2 + 2jh m v 2 l - 2j 2 l 2 , (45) 
* = !,..., LVJ, j = l,...,[hr, 2 \, 



and 



-Rij 



A+iji^ 2 + A+j' 2 , 



= 2ih m vik - 2i 2 k 2 + 2jh m v 2 l - 2j 2 Z, (46) 
i = 0,1, LVJ, j = l,...,[h m \. 

Obviously, these modes should satisfy a 3+1D Dirac 
equation. The above KG equations and the resultant 
equations for the profiles give clues for what the form 
of these should be. Each mode belonging to iff- or ty + 
has a z-dependent profile g™ 1 or y™ respectively. In turn, 
each chiral mode inside and ty + has a y-dependent 
profile /™ L R , and as we noted above, /™ L R = f™ L R . In 
using this expansion and performing dimensional reduc- 
tion, one expects a left-chiral mode of "J- with a par- 
ticular mass to not only attain a mass term with the 
corresponding right-chiral modes of <I>_ but also those 
of This is particularly important when considering 
that for a left-chiral mode in which has X^_ L = but 
A 2 7^ (ie. the mode is a zero mode with respect to 
the y-dependent wall) given that there is no correspond- 
ing right-chiral mode of the same eigenvalues and thus 
mass in However, there is in ty + and thus with such 
modes there is a single Dirac fermion of the given mass 
formed from the left-chiral mode in and the right- 
chiral mode of Hence, the correct ansatz is that the 
effective 3+1D mass Lagrangian should be of the form 



C 



h.c. 



For the massive modes with both A^r and \ 2 _ non- 



zero, one can deduce from the mass matrix in Eq. 47 that 



there exist two Dirac fermion modes of mass J A^ L + A?_ . 
After putting this ansatz into the 5D Dirac equation and 
doing some algebra, one can show that the equations 
yielded for the y and z-dependent profiles are exactly 
the same as those derived above from the Klein-Gordon 
equation. 

In this section we have shown that there exists a single 
chiral zero mode localized to the domain-wall intersec- 
tion when the Yukawa interactions in Eq. [21] In ad- 
dition to this single chiral zero mode, there is a single 
tower of Dirac modes which are zero modes with respect 
to one wall but not the other, and then for each given 
squared mass value for which the eigenvalues associated 
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with each wall are both non-zero there exist two Dirac 
modes. There will also be modes with one of \\ L R or \± 
being more than the maximum value for the localized KK 
modes and the other corresponding to a value associated 
with a localized mode; these modes can propogate along 
one wall and behave as 5D delocalized modes. Modes for 
which both the eigenvalues A± L R and Aj_ are more than 
the maximum values for the localized modes are com- 
pletely delocalized from both walls and can propogate 
through the entire 6D bulk. 



V. SCALAR LOCALIZATION 

Scalar localization can be similarly achieved via quartic 
coupling to the background scalar fields. Generally, in a 
model that is to be physically viable we are interested in 
localizing Higgs fields which have gauge charges. Thus 
for this section we will assume that our candidate scalar 
field is a complex scalar field The scalar potential for 
$ is then 



(48) 

Assuming that either $ has a vanishing vacuum expec- 
tation value (VEV) or one of much smaller magnitude 
than those that the background fields attain (as it would 
be in the case of an electroweak Higgs boson) , we can ig- 
nore the quartic self-coupling for <p in the determination 
of the profiles when we do a mode expansion. Hence, we 
can focus solely on the couplings of $ to the background 
fields and the mass term, and to calculate the profiles we 
must solve the 5+ ID KG equation 



[D+fi% + X^niVi + A *xiXi + X'S-mxiViXi 
+ Xt-mV2 + A * X 2X2 + A*^ X2 772X2] $ = 0. 



(49) 



Again assuming the same perpendicular solution for the 
background that we assumed in the previous section, we 
expand $ as a series of modes 



(50) 



Schrodinger equations 

[-- 
L dy 2 

+ ^<fvixiVi(y)xi(y)]p m (y) = X] n p m (y) (51a) 

d 2 

[ - ^2 + X ~fmVi(z) + >"5>x2xl{z) 

+ ^ mX2 -q 2 {z)x2{z)\q m {z) = X 2 n q m {z), (51b) 



2 <-Xi m vi{y) + ^xixi{y) 



reduces Eq. [49] to a relation between the masses of the 
KK modes m to the eigenvalues \} n , \ 2 m and the 5D bare 
mass /i$ 



m 2 = fi% + Xl n + Xl 



(52) 



Working in the non-dimensionalised coordinates y\ = 
y = ky, 2/2 = 2 = Iz, using the notation p} n = p m and 
Pm = Imi an d given the perpendicular solution given in 
Eq. |11[ the Schrodinger equations |51| can both be rewrit- 
ten in the form 



d 2 Pl 

dyf 



+ VHs(y~i)p l m (yi) = E lnPm(yi)i 



VhsiVi) = A 2 + (B 2 - A 2 - Aj) sech 2 (&) 
+B l {2A l + 1) sech (&) tanh (&), 

where the Aj and _B, are defined as 



(53) 



(2[(A$ Xi - - -) 2 + A| n . x . 



B, 



: 2Ai + r 



2> 



(54) 



and the eigenvalues of these potentials E % m are given in 
terms of the original eigenvalues X l m and Ai as 



E' 



Ai 



X$rn + A-i 



(55) 



Assuming A\ and A2 are positive 1 , these hyperbolic 
scarf potentials yield a discrete set of modes localized 
to the intersection of the domain walls. A localized <j) m 
mode must clearly have both of p m and q m decay to zero 
as y —¥ ±00 and z — > ±00 respectively. Since the respec- 
tive hyperbolic scarf potentials for the p m and q m yield 
\A{] and \A 2 ~\ localized functions respectively, there are 
\A{\ \Ai\ modes localized to the domain- wall intersec- 
tion. For each potential, the eigenvalues for the localized 
modes are known to be 



El 



2nA,, 



n 2 , 



(56) 



where the (j) m {x^) are 3+1D scalar modes satisfying the 
KG equation □<7) m (x Al ) = — m 2 </> m (a: A1 )and the p m {y) and 
q m (z) are the associated profiles along the y and z di- 
rections respectively. Subsitituting this expansion into 



Eq. 49 and then demanding that profiles satisfy the 



1 If one of A\ or Ai is negative, then there are no modes localized 
to the intersection and there will exist modes localized to one 
wall but delocalized from the other. If both are negative, then 
all modes are delocalized from both walls. 
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for n = 0, 1, [Ai\, and using Eqs. 52 55 and 56 we 
thus find that the squared masses of the localized 3+1D 
modes are 



/4 + \<s> m k + U V2 l - (^i - ifk z - (A 2 



-3?l\ 
(57) 

for £ = 0,1, L-4iJ and j = 0, 1, [A 2 \. 

There will also exist modes with sufficient average mo- 
menta transverse to one domain wall but not the other 
such that these modes are localized to one wall but not 
the other. These modes have a profile with one of the 
quantum numbers described above with respect to one 
wall but will have energies above that of the most en- 
ergetic localized mode of the other. These modes are 
essentially 5D particles. 

For modes with energies above all localized modes, 
their profiles along both directions are delocalized and 
they are thus 6D particles with can propogate along the 
full extent of the bulk. 

In any potentially phenomcnological model based on 
this type, we are usually interested in the lowest en- 
ergy localized mode, as this 4D scalar mode would cor- 
respond to Higgs particles in the effective field theory on 
the domain-wall intersection. In our model, this mode is 

The resultant (non- 



0, j = mode of Eq. 57 



the 

dimcnsionalised) profiles for this mode po{y) and qo(z) 
are simply given by 



Po{y) = C?,e~ Al log ' cosh to)! -2 - 81 arctan I tanh (5/2)] 

qo(5) = ^gg-^a log [cosh (z)]-2S 2 arctan [tanh (z/2)] 



(58) 



VI. CONCLUSION 

In this paper, we have generated a rare analytic solu- 
tion to a Z 2 x Z 2 -invariant scalar field theory with four 
real scalar fields in 5+ ID spacetime describing a pair of 
intersecting domain walls with internal structure. We 
found that with respect to the desirable perpendicular 
solution, there also existed a class of solutions describing 
kink-lump solutions which intersect at an angle between 
and 90 degrees as well as a solution where the walls are 
parallel and that these solutions were energy degenerate. 
We then argued that there exists a conserved topological 
charge related to the one-dimensional boundary of the 
y — z plane which differed between the intersecting solu- 
tions and the parallel solution, meaning that the inter- 
secting solutions cannot evolve to the parallel one despite 
this energy degeneracy. We also gave an argument as to 
why the perpendicular solution might be energetically fa- 
vorable to the solutions with intersection angle less than 
90 degrees in a nearby region of parameter space. 

In addition to finding this solution, we showed in the 
case that the two domain walls were perpendicular that 
fermions and scalars could be dynamically localized to 
the intersection of the two walls. We found that cou- 
pling a 6D fermionic field to one kink-lump pair with 
ordinary scalar Yukawa couplings and to the other with 



pseudoscalar Yukawa couplings allowed the fermionic sec- 
tor to be invariant under the full Z 2 x Z 2 symmetry and 
resulted in the localization of a 4D chiral zero mode on 
the intersection of the kinks, followed by a tower of lo- 
calized KK modes and 5D and 6D delocalized modes. 

Standard quartic couplings of a complex scalar field to 
the background scalar fields resulted in a tower of local- 
ized 4D localized scalar modes with the squared masses 
starting from some potentially non-zero value. This re- 
sult is similar to the result for scalar localization to a 
single domain wall in 5D. Furthermore, the squared mass 
of the lowest energy scalar mode can be negative, allow- 
ing the possibility of the localized scalar field inducing 
spontaneous symmetry breaking should that scalar field 
transform under a non-trivial gauge group representa- 
tion. 

Localization of gravity and gauge bosons is left to later 
work. Localization of gravity would involve searching 
for a similar solution to the 6D Einstein field equations 
as well as the Einstein-Klein-Gordon equations. These 
equations are highly non-linear and difficult to solve so it 
remains to be seen if an analytic solution could be found. 
It could also be the case that such a solution could have 
qualitative differences as in principle the two domain- wall 
branes should interact gravitationally, whereas in this flat 
space case the net interaction between the two kink-lump 
solutions was zero. 

For the localization of gauge bosons, we conjecture that 
the Dvali-Shifman mechanism works in 5+1D spacetime. 
This ultimately depends on whether or not non-Abelian 
gauge theories are confining in 5+1D. In previous work 
in 4+1D, particularly in the SU(5) model [T3], the Dvali- 
Shifman mechanism was facilitated by the addition of a 
scalar field transforming under the gauge group which at- 
tained a lump-like profile and induced symmetry break- 
ing in the interior of the domain wall. Hence, the back- 
ground solution of this model was of the kink-lump type 
discussed in this paper, and the additional fields xi and 
Xi which attained lump-like profiles in the 6D model are 
potentially well motivated. With the assignment of gauge 
representations, these fields can break a gauge group to 
different subgroups which then clash at the intersection, 
leading to further symmetry breaking and a different re- 
alization of the Clash-of-Symmetries mechanism [151 134] . 
We will discuss in detail this realization of the Clash- 
of-Symmetries mechanism and its applications to model 
building in a later paper. 
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